Production of the u meson in the pd ^He uj reaction at 1450 MeV and 1360 MeV 
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The production of u) mesons in the 



Hecj reaction has been studied at two energies near 



the kinematic threshold, Tp — 1450 MeV and Tp = 1360 MeV. The differential cross section was 
measured as a function of the u cm angle at both energies over the whole angular range. Whereas the 
results at 1360 MeV are consistent with isotropy, strong rises are observed near both the forward and 
backward directions at 1450 MeV. Calculations made using a two-step model with an intermediate 
pion fail to reproduce the shapes of the measured angular distributions and also underestimate the 
total cross sections. 



PACS numbers: 14.40.Cs, 25.40.Ve 



I. INTRODUCTION 



Extensive studies of oj production in the 7r~p — > noj 
reaction were carried out in the 1970s at the NIMROD 
synchrotron from the kinematic threshold up to uj cm 
momenta of = 260 MeV/c P, 0, Q. Though the dif- 
ferential cross sections were found to be isotropic, the 
data also showed a remarkable suppression of the w pro- 
duction amplitude near the threshold 0, Q . In the work 
of Binnie et al. where ir^p nrj', Tr~p n(f> and 
Tr~p —^nrj were studied with the same apparatus, no 
similar effects were found in the rj' and (j> cases and a 
threshold enhancement, rather than a suppression, was 
observed for the rj. First, a final state interaction (FSI) 
effect, where one of the pions from the u decay scat- 
ters off the recoiling neutron, was suggested as a possible 
explanation for the lu threshold suppression f]\ . This ex- 
planation was however tested and abandoned in Keyne et 
al. 0, where an alternative explanation were advanced 
in terms of a combination of s- and p-wave resonances. 



The data available for ui production in pd — > ^Hew 
are rather scarce. The reaction was first studied near 
its kinematic threshold at SATURNE 4]. The differen- 
tial cross section at Tp = 1450 MeV (p* = 280 MeV/c) 
measured in the forward and backward regions with the 
SPESIII spectrometer showed clear anisotropy in 6** , 
with strong peaking at extreme angles. This is in con- 
trast to the angular distributions observed in Tr^p —> noj, 
which remain fiat up to at least p* « 200 MeV/c 

Wurzinger et al. measured the pd ~^ ■^Hew cross sec- 
tion at 6** — 180° as a function of energy using the SPE- 
SIV spectrometer at SATURNE [6] and found a suppres- 
sion in the production amplitude near threshold similar 
to the one observed in the tt~p ^ noj reaction [J, H, [1] . 
These authors also described their data in terms of the 
FSI hypothesis. However, the correctness of this in- 
terpretation has been questioned for both reactions in 
Ref. 0, where effects associated with the uj width 
have been stressed. 

The production of heavy mesons in pd collisions has 



been studied theoretically in a two-step model d, [13, [HI, 
[12, [H, [13, [IE [III, which first involves the production 
of a light meson in the interaction between the incident 
and one of the target nucleons. In the second step, the 
light meson interacts with the other target nucleon to 
create the observed heavy meson. This procedure allows 
the large momentum transfer to be shared between the 
nucleons. The predictions of this model have been eval- 
uated for the specific pd ^^Keu reaction [iJL but only 
at the backward angles where data existed and no 
previous attempts have been made to calculate the full 
angular distribution. 

In this paper, the results of the measurements of uj 
production in the pd — ^ "^Hew reaction by the CEL- 
SIUS/ WAS A collaboration are presented at two differ- 
ent beam kinetic energies, Tp = 1450 MeV and Tp = 
1360 MeV. These correspond, respectively, to to cm mo- 
menta of p* = 280 MeV/c and p* = 144 MeV/c, i.e., 
excess energies of 63 MeV and 17 MeV. We have pre- 
viously reported data on the angular distribution of the 
UJ decay plane that show that the oj is produced largely 
unpolarized [l7| . which is in stark contrast to the al- 
most complete polarization of the (f) from the analogous 
pd '^He^ reaction [l3|. We now provide the corre- 
sponding results on the differential cross sections that 
have been measured over the full angular range. These 
new data constitute an important input to the intriguing 
question of how heavy mesons are produced in few-body 
collisions and, in particular, whether the reaction can be 
usefully viewed as a sequential two-step process. 

The outline of this paper is as follows. We start by in- 
troducing the WASA detector before going through the 
analysis in some detail. This involves principally the '^He 
identification, event selection, background, and normal- 
ization procedure. After describing the acceptance cor- 
rections, we show our results and compare them to calcu- 
lations that have been performed in terms of a two-step 
model. 



II. THE CELSIUS/ WAS A EXPERIMENT 

The measurements of the pd —> ■^Hew reaction were 
carried out at the The Svedberg Laboratory in Uppsala, 
Sweden, using the WASA detector which, until June 
2005, was an integral part of the CELSIUS storage ring. 
A deuterium pellet target [H, [2l| was used in combina- 
tion with a proton beam. The ^He were detected in the 
WASA forward detector (FD), which covers laboratory 
polar angles from 3° to 18° with respect to the beam di- 
rection. The forward detector consists of a sector-like 
window counter (FWC) for triggering, a proportional 
chamber (FPC) for precise angular information, a trig- 
ger hodoscope (FTH) for triggering and off-line particle 
identification, a range hodoscope (FRH) for energy mea- 
surements, particle identification and triggering, an inter- 
mediate hodoscope (FRI) for improved track reconstruc- 
tion and neutron detection and finally a veto hodoscope 
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FIG. 1: Side view of tlie CELSIUS/WASA detector setup. 
The Central Detector built around the interaction point (to 
the left) is surrounded by an iron yoke. The layers of the For- 
ward Detector are shown in the right-hand side. The CEL- 
SIUS beam pipe runs horizontally and the target pellets are 
injected through the vertical pipe. The abbreviations in the 
figure are explained in the text. 



(FVH) for triggering. All FD elements, except the FPC, 
are made of plastic scintillators. 

Mesons and their decay products are detected mainly 
in the central detector (CD), which consists of the Plas- 
tic Scintillating Barrel (PSB), the Mini Drift Cham- 
ber (MDC), and the central Scintillating Electromag- 
netic Calorimeter (SEC) that is made of Csl crystals. 
Charged particles, such as pions from uj decay, are dis- 
tinguished from neutral particles by their signals in the 
PSB, which also provides angular determination from 24° 
to 159°. The momenta of charged particles are estimated 
by tracking in a magnetic field from the SuperConduct- 
ing Solenoid (SCS) using information from the MDC. The 
SEC measures angles and energies of photons from me- 
son decays in the polar angle range from 20° to 169°. 
A schematic overview of the CELSIUS/WASA detector 
setup is shown in Fig. [T] 

A special trigger was developed to select "^He events, 
which are characterized by one high-energy-deposit hit 
in the FWC and one hit that overlaps in the azimuthal 
angle in one of the consecutive detectors. 

III. ANALYSIS 
A. ^He identification 

The '^He ions are identified in the FD by the AE — E 
method. For this the scintillation light output in the de- 
tector layer where the particle stops is compared to that 
in the preceding layer. The particles then show up in 
different bands, as illustrated in Fig. [2l where the up- 
permost one represents the '^He. For these candidates 
the scintillation light output is converted to energy de- 
posit, taking into account the light quenching [2^. This 
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FIG. 2: The scintillation light output in the first layer of the 
FRH versus the light output in the second layer for parti- 
cles stopping in the second layer. The uppermost band corre- 
sponds to ^He ions, below which is a band of stopping protons. 
The intense spot at the bottom left corresponds to minimum 
ionizing particles. 



FIG. 3: The acceptance for the pd — > ''He a;, uj — > tt^tt^tt'' 
reaction as a function of cos 0^ at 1360 MeV (solid line) and 
1450 MeV (broken line) after applying the cuts explained in 
the text. The dip in the 1450 MeV distribution corresponds 
to the ^He that stop between the first and second layer of the 
FRH. 



allows the corresponding range to be estimated, which is 
added to the amount of material traversed by the particle 
in the preceding detector layers to give the total range. 
Energy-range tables then allow the initial kinetic energy 
to be estimated. The expected light output in all the 
detector layers traversed is calculated and compared to 
the measured one by evaluating the x^, according to 

N 

X^=^(ALr-AL^)Va?, (1) 

1=1 

where i is the layer number, N the number of layers tra- 
versed, AL™ the measured light output, AL1 the calcu- 
lated light output, and Oi the light output uncertainty 
in a given layer. The particle is then considered to be 
properly identified if the calculated does not exceed 
Xmax = 6.0. For more details, see Ref. [23j . 

B. Event selection 

In this work, we focus primarily on the uj — > 
tt'^tt^tt" decay channel because the large branching ra- 
tio (BR=89.1%) gives the high statistics required for the 
extraction of angular distributions. However, a parallel 
analysis of the uj 7r°7 decay channel (BR=8.7%) is a 
valuable tool for checking that cut efficiencies and other 
experimental biases are under control. 

1. pd — > ^Heuj, ui — > ■k^-k'^-k~ events 

To select u> — >• tt^tt+tt" events, we require the ^He ion 
to be measured in the FD with well defined energy and 



angle, which means that the "^He has to stop within the 
FRH. The geometrical acceptance of the FD is 95% at 
1450 MeV and 78% at 1360 MeV. The main event loss is 
due to the '^He emitted at small angles that escape down 
the beam pipe. The detection efhciency of the FD is fur- 
ther reduced by nuclear interactions in the detector ma- 
terial. The efficiency of detecting a ^He from pd — > ^He w 
in the FD, as deduced from Monte Carlo simulations in- 
cluding geometry and detector material responses, is 61% 
and 54% at 1450 and 1360 MeV, respectively. 

In addition to the ^He selection, at least two pho- 
tons are required in the SEC. Furthermore, one photon 
pair must have an invariant mass that does not differ 
from that of the 7r° by more than 45 MeV/c^, and the 
missing mass of the '^He system must be larger than 
250 MeV/c^, i.e., twice the pion mass, after taking the 
resolution into account. Finally, two or more hits are 
needed in the PS and the total energy deposit in the 
SEC must not exceed 900 MeV. This rejects candidates 
with photon signals from accidental events. 

These constraints give an overall acceptance of 35% 
at 1450 MeV and 34% at 1360 MeV. The differential 
acceptance is shown in Fig. [3] for the two energies as a 
function of cos0* . 

2. pd ^Heco, LU -k'^j events 

All final state particles, i.e., the '^He and three photons, 
can be measured with high acceptance in the uj — > tt'^j 
case. Each event is fully reconstructed, which ensures 
that the kinematical constraints are fulfilled. Though 
this leads to a cleaner sample than the three-pion chan- 
nel, the low branching ratio (BR=8.7%) yields statistics 
that are insufficient to provide angular distributions. 
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Having identified a ^He and three photons, we de- 
mand that one photon pair has an invariant mass close 
to that of the tt", |/Af(77) - m^n| < 45 MeV/c^, and 
that the invariant mass of all three photons is larger than 
600 MeV/c^. The missing-mass-squared of the ■^He777 
system must not exceed (100)^ MeV^/c'*. We also re- 
quire that the difference between the directions of the 
missing momentum of the '^He and that of the 87 sys- 
tem be smaller than 20°. As in the three-pion case, a 
limit is placed on the total energy deposit in the SEC, 
which may not exceed 1200 MeV. Finally, a coplanarity 
cut is applied: 160° < |(/>iab(^He) - 0iab(37)| < 200°. 
These constraints reject very effectively the contribution 
from accidentals and events originating far from the tar- 
get region. They give acceptances of 19% and 18% at 
1450 MeV and 1360 MeV, respectively. 



from the pd ^HcTr^Tr" reaction, where one of the 
photons escapes detection. Assuming phase space pro- 
duction, about 1.8% of the 27r° events survive the cuts 
optimized for 7r°7 selection at 1450 MeV and 1.4% at 
1360 MeV. Although these are tiny fractions, the two- 
pion production cross section is large and the uj tt'^j 
branching ratio small. Hence it is expected that the num- 
bers of events surviving the cuts should be of the same 
order of magnitude for the signal and background. 

Interactions between the beam halo and the beam pipe, 
as well as chance coincidences, may also contribute to the 
background. However, we show in Sec. IIVI that they only 
do so to a very limited extent. 

D. Normalization 



C. Backgrounds 

In the three-pion case, the principal background is di- 
rect pion production, i.e., pd — > '^He7r''7r+7r^, which has 
exactly the same signature as uj production. The accep- 
tance for this background, assuming phase space produc- 
tion, is 31% at 1450 MeV and 35% at 1360 MeV. 

Both Lo and direct Stt events are produced mainly 
through the interaction of the beam with the pellet tar- 
get. However, such reactions can also occur through in- 
teractions of the beam with the rest gas, the beam halo 
with the rest gas, and the beam halo with the beam pipe. 
Rest gas isproduced when the pellets are vaporized by 
the beam (2J] or when the pellets are not properly cap- 
tured in the pellet beam dump. In the present work, we 
found that around 30% of the data was not produced 
from the region of the pellet target. This is in line with 
the figures obtained from a recent measurement of the 
pd ^He77 reaction at Tp = 893 MeV The vapor- 
ization of the pellet is expected to be higher at lower en- 
ergies |24i] , but the amount of rest gas also depends upon 
other time-dependent experimental conditions, such as 
the beam-target overlap, the efficiency of the vacuum 
pumps, and the performance of the pellet dump. Since 
the energy reconstruction procedure assumes that all par- 
ticles come from a well defined interaction point, the 
measured energy and momentum of particles produced 
far from the pellet target will have larger uncertainties 
and thus degrade the missing-mass resolution. However, 
the width of the ui peak agrees well with Monte Carlo 
simulations for a well defined interaction point, which 
means that most events come either from the target or 
from the close vicinity. The latter can be treated as good 
events that do not significantly affect the quality of the 
data sample. 

Monte Carlo simulations show that other channels, 
such as pd - 
^He TrOTrOTT+TT 



^HcTT+TT , pd '^He7r°7r° and pd 



give negligible contributions to the back- 
ground and will be ignored. 

The main background channel in the ir^j sample comes 



The luminosity depends on the beam intensity and the 
pellet rate, which are monitored and recorded during the 
data taking. However, the luminosity also depends on 
the overlap between the beam and target, which cannot 
be measured directly. It is also not constant during a 
run or even over a cycle of data taking. Furthermore, 
as previously mentioned, as many as 30% of the events 
may result from interactions away from the target. This 
means that the luminosity must be evaluated by normal- 
izing the number of events of some subsidiary reaction 
to a well known cross section. Quasi-elastic pp scatter- 
ing might be a good choice because the cross section is 
well known and large statistics are available. However, 
these data were collected with a different trigger from 
the '^He events and the relative trigger efficiency is not 
known with sufficient precision. It is therefore preferable 
to use a reaction with a "^He in the final state, where the 
trigger was the same as for pd ^ ^Heuj. 

The differential cross section for the pd ^Herj reac- 
tion was measured at a few very backward angles with the 
SPESIV spectrometer at 1450, 1350, and 1250 MeV 0]. 
It was also measured over a larger angular range at 
1450 MeV with the SPESIII spectrometer The two 
data sets agree well in the region of overlap. WASA has 
collected clean pd "^Her;, rj samples at 1450 

and 1360 MeV. In addition, we analyzed the pd ^Yiorj 
reaction via the 77 7r°7r~''7r~ decay channel, which pro- 
vided a cross check of the cut efficiency. There is a good 
agreement between the number of rj mesons in the 27 
and 37r samples This shows that the effects of the 

cuts are well understood and that the yield of the events 
coming from outside the target is about the same in the 
pd '^Hc77, rj ^ case as when the cuts are optimized 
for selecting a "^He tt^tt+tt" final state (including tt'^tt+tt^ 
production via rj and uj). 

The integrated luminosity L is calculated from: 
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where AA'' is the acceptance-corrected number of 77 can- 
didates in the angular region Afl = 27rA(cos0*), and 
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FIG. 4: Differential cross section data for the pd '^Her; 
process at 1450 MeV, which has been used for the evaluation 
of the luminosity. The filled and unfilled triangles show data 
from SPESIV 26]] and SPESIII 0, respectively. The line is a 
constant fitted to both sets of data points. The filled circles 
represent the WASA data after normalization. 
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FIG. 5: Differential cross section data for the pd — > '^He?; re- 
action used for the evaluation of the luminosity at 1360 MeV. 
The unfilled triangles are data from SPESIV ^26|] taken at 
Tp — 1350 MeV; the point at cos 6^ = —0.96 is used for 
normalization. The filled circles are WASA data taken at 
1360 MeV and the line is the result of a linear fit to these. 



IV. RESULTS 



da/dO, is the differential cross section at this angle. 

At 1450 MeV, there is a region -0.99 < cos 6** < -0.6 
where two WASA points overlap with a set of SPE- 
SIII and SPESIV points (see Fig. [4|). In this region, 
the angular distribution is approximately isotropic. A 
fit to the combined SPESIII and SPESIV points to a 
constant gives a x^/ndf = 1.8. The resulting cross sec- 
tion, da/dO. — 2.53 ± 0.19 nb/sr, where the statistical, 
systematic and normalization uncertainties of each point 
have been taken into account. A fit to the two WASA 
acceptance-corrected numbers of counts to a constant has 
an uncertainty of 9% and this contributes to the 12% 
overall uncertainty in the resultant luminosity. 

The poorer pd — > ^He rj database at 1360 MeV makes 
the situation more difficult at this energy and we have to 
rely on the few SPESIV points taken at Tp = 1350 MeV 
in the near-backward region [2^1 ■ At 1250 MeV and 
lower, this data set extends slightly further in angle and 
from this it is clear that, in addition to the very sharp 
backward peak, the angular distributions show a mini- 
mum for COS0* « —0.96 before rising less steeply at more 
forward angles. Though the four WASA points in the 
backward hemisphere shown in Fig. O are consistent with 
such a small rise, they do not overlap with the SPESIV 
data so that a linear fit has been made to extrapolate 
to the average of the two points in the cos 6*,* = —0.96 
region. It could be checked from the SPESIV data at 
1450 MeV that the energy dependence in this region is 
negligible compared to a total normalization uncertainty 
of 29%, of which 27% arises from that in the linear ex- 
trapolation. 



A. Tp = 1450 MeV 

The u} mesons are identified from the peak in the '^He 
missing-mass distribution. The number of uj events is ex- 
tracted by subtracting the background in two ways. The 
first is by fitting the background to a phase space Monte 
Carlo simulation of pd ^He tt^tt+tt" data. The sec- 
ond is by fitting the data as a Gaussian peak sitting on 
a polynomial background. The number of background 
events can then be estimated by integrating the polyno- 
mial. The central value and the width of the Gaussian 
peak obtained in the fit provide a quality control. 

The difference in the number of lu events extracted 
by subtracting the background in these two ways gives 
the largest contribution to the systematic uncertainty. 
Applying the constraints defined in Sec. IIIIB II gives a 
sample with a '^He missing mass distribution shown in 
Fig. [HI After subtracting the background, we are left with 
9900 ± 700 candidates for the u; — > tt^tt'^tt^ channel. 

To verify that the systematics are under control, we 
made a cross check using the uj — > n'^j channel. On 
the basis of the evaluated acceptances and the known 
branching ratios, both summarized in Table HI we expect 
520 ± 40 events from oj — > 7r'^7 in the sample. Selecting 
the events according to method described in Sec. IIII B 21 
we get a number of 420 ± 50, which is fairly consistent. 

The good match between the background continuum 
and the simulated pd ■^HeTr^Tr+Tr" in Fig. [HI the good 
agreement between the w width in Monte Carlo and data, 
and the consistency between the results for the 37r and 
tt'^j channels gives us confidence that the contribution 
from background not coming from reproducible physical 
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MM( ^He) in MeV/c^ 

FIG. 6: Upper panel: Points are the distribution of missing 
mass of the ^He, MM(^He), for all events at 1450 MeV fulfill- 
ing the constraints optimized for pd ^fieco, uj — » tt^tt^tt" 
selection. The histogram shows the corresponding phase 
space simulations of the pd — > ■'He7r''7r"'"7r~ channel. Lower 
panel: MM(''He) distribution for events fulfilling cuts opti- 
mized for selection of the pd "^Hetj, uj — > iv'^'y channel. The 
histogram is the phase space simulation of pd — » ^HeTr'^Tr'^. 
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TABLE I: Summary of the numbers of lj candidates found at 
different energies for the two decay channels. 



Channel 


Tp 
(MeV) 


Acceptance 
(%) 


BR (%) Number of 
(%) candidates 


u> — > 7r'''7r~7r" 


1450 


35 


89.1 


9900±700 


U! — > 7r"7 


1450 


19 


8.7 


420±50 


+ - u 

ClI — > TT TT TT 


1360 


34 


89.1 


1800±200 




1360 


18 


8.7 


80±20 



reactions may be safely neglected. 

In order to extract the differential cross section for the 
Stt sample, we divided the data into cos 0* bins and plot- 
ted separate "^He missing-mass distributions, some ex- 
amples of which are shown in Fig. [71 The number of uj 
events was then obtained by subtracting the background, 
which was estimated in the two ways described for the to- 
tal sample. The systematic uncertainty from background 
subtraction varies between 5% and 20%. Adding the uj 
candidates in each bin gives 9090, which is in good agree- 
ment with the value quoted in Tabic [H The numbers 
were corrected for acceptances that were derived from 
the phase space Monte Carlo simulations. Since the an- 
gular distribution is highly anisotropic, a fitted second 
degree polynomial was used as input. This reduced the 
total acceptance from 35% to 33%, though the effect for 
individual bins was negligible. This shows that the model 
dependence of the acceptance on angle is very weak. Fi- 
nally we converted the numbers of uj mesons to cross 
sections by using the integrated luminosity discussed in 
SeciniDl 



FIG. 7: Data taken in four angular regions at 1450 MeV, with 
cuts optimized for the pd — > ^Heo;, u! — > tt^tt^tt" selection. 
Phase space simulations of pd ^HeTr^Tr+Tr" production are 
also shown, as well as a fitted Gaussian peak on a polynomial 
background. 



The resulting pd — > ■^Hew angular distribution is 
shown in Fig. [S] along with the values obtained at SPE- 
SIII 0. The two data sets agree in the backward direc- 
tion, which lends confidence to our normalization pro- 
cedure. However, although the WASA data do show a 
slight increase towards the forward direction, the sharp 
peak observed in Ref. is not confirmed. 

Fitting the angular distribution with Legendre poly- 
nomials results in an integrated cross section of ctot = 
83.6±1.5±2.2 nb, where the first uncertainty is statistical 
and the second systematic. In addition, however, there 
is an uncertainty of 12% that comes from the luminosity 
determination. 



B. Tp = 1360 MeV 

The analysis is more difficult at 1360 MeV than at the 
higher energy. The finite uj width (F = 8.44 MeV/c^) 
leads to an asymmetric peak close to threshold since the 
high mass tail of the uj meson is then partially suppressed. 
Secondly, the background subtraction is more difficult 
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FIG. 8: pd ^^Ueto differential cross section at 1450 MeV. 
The error bars on the WASA data (filled circles) show the 
statistical uncertainties while the grey histogram shows the 
systematical uncertainties due to background subtraction and 
acceptance correction. In addition there is an overall normal- 
ization uncertainty of 12%. The SPESIII results are shown 
by the open triangles. 



since the multipion continuum ends under the lu peak 
(see Fig. Thirdly, the signal-to-background ratio is 
small. To add to these difhculties, fewer data were taken 
at 1360 MeV than at 1450 MeV due to beam-time con- 
straints. Higher systematic and statistical uncertainties 
are therefore to be expected at this energy. 

Selecting data according to the methodology of 
Sec. nil B II gives the data sample shown in the upper 
panel of Fig. [5] which, after background subtraction, 
leaves 1800 ± 200 u; events. Using the information given 
in Table HI we then estimate the number of w ^ 7r°7 
events should be 90 ± 10. This is to be compared to 
80 ± 20 events obtained by subtracting the background 
in the lower panel of Fig. [9l The numbers from the two 
decay channels are thus quite consistent. 

The extraction of the angular distribution was carried 
out in the same way as at 1450 MeV. The sum of the 
events in each angular bin (1600) agrees well with the 
1800 found from Fig. [51 The angular distribution shown 
in Fig. [TO] is consistent with isotropy and, within the to- 
tal uncertainties, agrees with the point at 1360 MeV ob- 
tained by interpolating SPESIV results Q. 

It was suggested in Ref. ^ that rescattering of de- 
cay pions off the "^He nucleus might explain the ob- 
served threshold dip in the production amplitude. Their 
Monte Carlo simulations suggested that this would lead 
to a difference in the measured cross sections for the 
uj — > tt'^tt+tt^ and uj 7r*^7 modes of around 10% at 
1360 MeV. Unfortunately, the uncertainties in the num- 
ber of ll! events from the two different channels obtained 
in the present experiment are too large to test this pre- 
diction. 



% 200 




MM( ^e) in MeV/c^ 

FIG. 9: MM(^He) data at 1360 MeV for both the Stt and 7r°7 
selection with experiment and simulation, as described for the 
1450 MeV data of Fig. H 
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FIG. 10; Angular distributions for the pd —+^¥LeLo reaction at 
1360 MeV, as described in Fig. [8] for the 1450 MeV results. 
In addition to the systematic uncertainties shown by the grey 
histogram, there is an uncertainty from normalisation of 29%. 
The triangle shows the differential cross section obtained by 
interpolating the SPESIV data [(|. The solid line represents 
the results of the model calculations described in Sec. [Vl 



—0.55 in Fig. [TOl which have especially large system- 
atic as well as statistical uncertainties, the angular dis- 
tribution is fairly flat and is in good agreement with 
the SPESIV point 6]. Excluding these points from a 
linear fit, the total cross section is determined to be 
CTtot = 84.6 ± 4.0 ± 4.8 nb. To this must be added the 
luminosity uncertainty of 29%. 



V. MODEL CALCULATIONS 

The minimum momentum transfer for the pd — > "^He w 
reaction at 1360 MeV is 1110 MeV/c and this only 
drops to ~ 935 MeV/c at 1450 MeV. A direct produc- 
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tion mechanism with such a large momentum transfer is 
expected to give vanishingly small cross sections. In or- 
der to share this large momentum between the nucleons, 
Kilian and Nann [28] suggested a two-step model, where 
a virtual pion beam is produced via a, pp dn^ reaction 
on the target proton, followed by the production of the 
observed meson X through the 7r+n — > pX reaction on 
the target neutron, as illustrated in Fig. [TlJ There is an 
analogous contribution with an intermediate vr". 




FIG. 11: Two-step model for the pd — > HeX reaction where 
a virtual pion beam, created through pp dn^ , produces 
the observed meson X through a second yr^n pX reaction. 
There is a similar graph with a neutral intermediate pion. 



The quantum-mechanical implementation of Fig. [TT] 
describes well the total cross section for the production 
of the 77-meson in pd collisions near-threshold [1, [ill, 13 ■ 
The ry angular distributions have been measured at sev- 
eral excess energies above about 20 MeV [l^. Although 
it had been claimed their shapes could be reproduced in 
a simplified two-step model llal. the most realistic im- 
plementation of this approach [l2| fails to reproduce the 
data. 

Theoretical studies of the pd — > '^Heo; cross section 
using a two-step model have been undertaken only close 
to threshold [lO, [3| and these underestimate the exper- 
imental data 6j. No predictions have yet been made for 
the LO angular distribution. 

Following closely the implementation of the two-step 
approach described in Ref. |_ll|, we evaluate the TriV 
LuN sub-process using the Giessen model |3fl] . This is an 
effective Lagrangian approach that takes seven coupled 
channels in account, viz., jN, ttN, 2-kN , rjN, ojN, KA, 
and KT,, for the simultaneous analysis of all the data up 
to 2 GeV in terms of 11 isospin-i resonances. It was 
shown in Ref. [sot that s-waves alone were insufficient 
to describe the ttN — > cuN data even in the very near- 
threshold region. The full data set required partial waves 
up to £ = 3 and their full amplitude analysis has been 
used as the basis for the construction of the vrA^ ^ uiN 
i- matrix that is the required input for the pd — > ^Hew 
estimation. The pp — > di:^ input was taken from the 
parameterized t-matrix of the SAID group [3]| . 

The results of the calculations performed in the plane 
wave approach to the two-step model, are shown in 
Figs. [TO] and [TO] at beam energies of 1360 MeV and 
1450 MeV, respectively, along with the available data. 
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FIG. 12: Comparison of the WASA data on the pd ^Reu 
reaction at 1450 MeV (filled circles) with calculations based 
upon a two-step approach (solid line). The meaning of the 
error bars and histogram is as in Fig. [S] 



Although the 1450 MeV predictions are reasonably close 
to the data in the central region, the shape is quite differ- 
ent to that of the experiment, giving forward and back- 
ward dips rather than peaks. There are therefore dis- 
agreements of up to an order of magnitude at the ex- 
treme angles. On the other hand, at 1360 MeV the two- 
step model predicts a fairly flat distribution with only a 
slight forward enhancement. However, it underestimates 
the cross section by a factor of about five in the forward 
region and by even more in the backward direction, which 
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is a similar discrepancy to that reported in Ref. 

It should be stressed that the pd ^Heui angular dis- 
tribution is predicted to be anisotropic at 1360 MeV even 
if only s-wave are retained for the TriV —^ luN vertex. This 
anisotropy persists in calculations carried out even closer 
to the threshold. It is perhaps interesting to note in this 
respect that experimental data on the pd — > '^He 77 reac- 
tion very near threshold show significant anisotropy [33 | . 

We compare in Fig.[T3]the total cross section calculated 
using the two-step model with the WASA data. Since this 
theoretical approach underestimates the differential cross 
sections, it also underpredicts the total cross sections. 

Why does the two-step model fail for uj production? If 
the intermediate pion is taken on-shell, as in the Kilian 
and Nann classical model [1^, the relative velocities of 
the deuteron from the pp — > dir'^ and the proton from 
the TT+n —>■ puj subprocesses are very large at the extreme 
angles at 1450 MeV. In contrast, at 1360 MeV, the change 
in the velocity matching over the full angular region is 
much less. The cross sections calculated using the two- 
step model at 1360 MeV thus describe the behavior of the 
experimental angular distribution but not its magnitude. 

The poor velocity matching implies that the interme- 
diate pion is a long way off its mass shell for both beam 
energies, unlike the case of the pd '^Herj reaction at 
threshold. The predictions might be improved by using 
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FIG. 13: The total cross section of pd '^Hew reaction as 
a function of the beam kinetic energy. The filled circles are 
the two WASA points, where the error bars include the nor- 
malization uncertainties. The solid line shows the estimation 
within the two-step model presented in this work. 



off-shell t-matrices for both pp dn^ and n^n puj. 
However, it is clear that to reproduce the data at the 
extreme angles at 1450 MeV requires contributions from 
other diagram(s) since poor velocity matching is a kine- 
matical effect that depends only upon the assumption 
of an intermediate pion. A similar conclusion might be 
drawn from the failure of the same model to describe well 
the pd 77 angular distributions away from thresh- 

old tli. 

VI. CONCLUSIONS 

The differential cross section for the pd ^Hecu re- 
action has been measured over the full angular range 
at 1450 MeV. The data show clear anisotropy, with 
strong rises in both the backward and forward directions. 
For the backward-going uj the agreement with SPESIII 
data [5*1 is good, but the sharp forward peaking claimed 
at SPESIII is not confirmed by our data. The two-step 
model fails to describe the angular distribution, giving a 



convex rather than the concave shape found in our data 
The angular distribution at 1360 MeV is consistent 
with isotropy and, within the experimental uncertainties, 
the differential cross section for the backward-going u 
agrees with the SPESIV point at cos 6** = -1 [1]. If 
the assumptions in Refs. [3, [1] that the SPESIV data 
are incorrectly interpreted are right, their point shown in 
Fig. 1 101 would be higher and the agreement with our data 
would become worse. 

The large angular coverage of the WASA detector al- 
lowed values of the total cross sections at both energies 
to be extracted for the first time. At 1450 MeV we find 
CTtot = 83.6 ± 1.5 ± 2.2 nb, with an additional uncertainty 
from the normalization of 12%. The corresponding num- 
ber at 1360 MeV is atot = 84.6 ± 4.0 ± 7.3 nb, where the 
additional normalization uncertainty is 29%. 

The two-step model underestimates the total cross 
section data by a factor of about two at 1450 MeV and 
by a factor of ten at 1360 MeV. This is probably due to 
the velocity matching between the proton and deuteron 
produced in the two subprocesses being poor for some uj 
angles. Further theoretical work is therefore needed to 
describe both w and r/ production in pd collisions. 
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